plot([w(1,1:2)],[w(2,1:2)])
plot([w(1,1:2:3)],[w(2,1:2:3)])

 t=0:0.01:2*pi;
 x=100*cos(t)+1;
 y=100*sin(t);
 plot(x,y,'g')
 axis equal